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Abstract 

Exact solutions of two-dimensional hydrodynamics equations for 
the symmetric configurations of two and four vortices in the presence 
of an arbitrary flow with a singular point are found. The solutions de- 
scribe the dynamics of the dipole toroidal vortex in accretion and wind 
flows in active galactic nuclei. It is shown that the toroidal vortices 
in a converging (accretion) flow, being compressed along the large ra- 
dius, are ejected with acceleration along the axis of symmetry of the 
nucleus, forming the components of two-sided jet. The increment of 
velocities of the vortices is determined by the monopole component 
of the flow only. The dipole component of the flow determines the 
asymmetry of ejections in the case of an asymmetric flow. 

PACS: 47.32.-y; 47.10.Df; 98.54.Cm 
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1. Introduction 



A large number of works (see, for example, the monograph [1]) is devoted 
to the origin of jets. In the most of them the decisive role is played by a 
strong magnetic field [2-4], or "external", or arising from the development of 
instabilities in the plasma of the accretion disk. This field serves as a guide 
for the movement of particles under the action of electromagnetic, centrifugal 
and gravitational forces, allowing them to move against the gravity and carry 
away the angular momentum that is necessary for the effective accretion 
process, which is responsible for the activity of the nucleus (see discussion 
and references in reviews [5-7]). At the same time, the very possibility of the 
existence of strong magnetic fields in the accretion disks around the black 
holes is not entirely clear. In this connection, the models of continuous flow 
without magnetic field are also considered, including those, in which the 
structure of jets resemble the hydrodynamic tornadoes [8]. 

The observations, however, show that at small distances from the nu- 
cleus (parsec scales for active galactic nuclei (AGN)) the emergence of the 
individual (including the superluminal) components of the radio jets are ob- 
served [9]. In the model we discuss here [10-12] some ejections are derived 
from the kinematics of the interaction of vortices and the exposure of the 
magnetic field is not required. An important role in this process has the flow 
that can affect the velocity of ejection. Therewith, greater velocities of the 
ejected components are attained in a converging (accretion) flow. 

We regard a system of toroidal vortices [10], surrounding the central part 
of AGN, the outermost of which is observed as a "obscuring torus". The 
vortex motion arises in the torus due to twisting by the wind and radiation. 
Due to the flow symmetry the movement in the torus possesses the dipole 
character (Fig. 1 in [10]). In its simplest form, this movement can be repre- 
sented as the motion of two opposite rotating vortex rings in a radial flow. 
As is known, the dynamics of a vortex ring can be described as the movement 
of a pair of point vortices that arise in the cross-section of the ring (torus) by 
the plane of symmetry. In our case it is a symmetrical system of two or four 
vortices (or two vortex pairs) in an arbitrary flow with a point singularity. 
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2. Setting of a problem 



We consider the motion of a system of point vortices on a background 
flow caused by the stationary singular point, which is placed at the origin. 
Following [13, 14], the stream function can be represented as the sum of two 
constituents: its regular part ip reg which describes the background flow and 
the singular part ip S ing, which describes the point vortices. 

The complex potential of the flow can be represented as a series 

w reg = Co In z H h — + ... (z = x + iy) (1) 

z z 

When the flow is symmetric about the axis Ox, all the coefficients C k in 
(1) are real. Expressing the complex potential through the usual velocity 
potential tp and the stream function ip according to w — tp + iip and passing 
on to the polar coordinates z — r • exp we obtain the stream function in 
the form 

^ req = C 9-— sin0- ^|sin20- .... (2) 

The radial and azimuthal velocity components of the background flow are 
determined by the conditions 

_ 1 dj) reg _ dj) reg 
r 89 ' dr 
Combining the monopole at the origin (source or sink), dipole, quadrupole, ... 
with intensity C ,Ci,C 2 , ... we can obtain a given distribution of the radial 
velocity component at the circle | z \— R which describes the considered flow. 
The system of vortices and the flow has Hamiltonian form. 

3. Symmetric motion of two vortex pairs in the flow with 
singularity of the type "source + quadrupole + ..." 



Now we consider the dynamics of a system of two vortex pairs in the 
stream flow generated by a fixed singular point provided the existence of two 
axes of symmetry (fig.l). This case can be interpreted as the motion of a 
point vortex in the right angle (which sides are impermeable "walls") at the 
apex of which is placed the referred hydrodynamic singularity. Note that 
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Figure 1: Scheme of the movement of four vortices (two vortex pairs) in 
a symmetric flow, containing the sink and the quadrupole for Cq = —1, 
d = 0, C 2 = l, T = -4vr 

the solution of the problem for a symmetric system of four vortices in the 
absence of the background flow was found in the classical Grobli's work (see 
[15, 16]). In a purely radial flow it admits a Hamiltonian formulation [11] 
and the exact solution of the dynamic problem [11, 12]. 

The stream function for the discussed flow is represented in the form 

VVe 5 = G)0-^f Sin20-^sin40-.... (3) 

The vortex components of the right pair are arranged symmetrically about 
the axis Ox in the points (x; y) and (x;-y), and the intrinsic (due to the 
interaction of vortices and not related to the presence of the background 
flow) velocity of the vortex in the 1-st quadrant is 

V*. r f 1 V . x 1 \ 

smg 47T [ y x 2 + y 2 ' x 2 + y 2 x J 

In the case of only one vortex pair without a background flow that corre- 
sponds to the known expression 

f s = (r/4ny;0), (4) 
where T is the intensity of the vortex, x and y are its abscissa and ordinate. 
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The dynamics of two pairs of vortices in the flow is described by the 
equations 



Air \y x z + y z J ay \x 2 + y z x) ox 

and the Hamiltonian of the system of 4 vortices is reduced to the form 

r 3j]J y *2<xy 

H = — In = + C arctan C2-^ ^2--- ( 6 ) 

4tt y/x 2 + y 2 x (x 2 + y 2 ) 

Accordingly the equations (5) can be represented as 

dH dH 



x = y 



dy ' dx 

Equation H = E = const determines the trajectory of vortices. Moreover, if 
the motion is unbounded, then the vortex pair comes from infinity along one 
axis {Ox), exchange their components and goes away to infinity along the 
other axis (Oy): 

f47r£l f A<k{E - C f ) 1 
yoo = exp |-p-| , x -> +°°; ^oo = exp <^ \ , y ->■ +oo; 

Co ^ 0, d ^ 0, C 2 ^ 0. 

For the asymptotic values and t/oo coordinates of the vortex moving in 
the first quadrant, we have the correlation 

yoo = Xooexp (27r 2 C /r) , (7) 

which coincides with the result for the purely radial flow u > 12 . This means 
that the ratio yoo/xoo of the limit distances between elements of the vortex 
pairs at infinity and, correspondingly, the ratio of their velocities, are only 
determined by the intensity of the C of the source (sink) at the origin and 
do not depend on the other multipole components. 

Since the velocity of the translational motion of a vortex pair is inversely 
proportional to the distance between them (cf. (4)), then we obtain the ratio 
of the asymptotic values of the velocities of the pairs, coming from infinity 
V_ = r/(47ry 0O ) and going to infinity V + = r/(47rx 00 ) 
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Figure 2: Phase portrait for the vortex of the first quadrant, corresponding to 
the movement of the four vortices in the accretion flow. Cq = — 1, C\ = 0, 
G% = \, T = —Att. Thin and bold lines correspond to the trajectories of the 
vortex and the separatrix accordingly. The asymptotic values of coordinates 
(Xoo and y^) correspond to half of the distance between the components of 
the pairs of vortices at the inlet and outlet of the system. 

V-/V+ = Wz/oo = exp (-27T 2 C /r) , (8) 

that coincides with the ratio obtained in [11, 12]. Obviously, the vortex 
pairs go away from the source at a slower velocity, and from the sink - 
corresponding to accretion - with a greater velocity than they come in. Thus, 
the obtained solutions confirm the main result of the work [12] about the 
acceleration of ejections by the radial accretion flow, and expand it on the 
general case of a symmetric accretion flow. As we will show below the dipole 
component of the flow may be responsible for the asymmetry of ejections. 

4. Motion of two vortex pairs in the flow from dipole which axis is 
the axis of symmetry of the flow 

In the case of motion of two vortex pairs in the flow from a dipole with the 
axis Oy with intensity C\ (that corresponds to replacement G\ — > iC\ in 
(1)) the equations of dynamics of the vortices can be expressed through the 
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coordinates of the two vortices in the right half-plane in the form 
dH dH dH dH 



a?/i ox i ay 2 0%2 

H = T 2 - T 2 + T~ ln ~ x 2) 2 + (2/1 - V2) 2 ) - (9) 

r / 2 2\ r r 

— In \{xi + x 2 ) + (t/i - 7/2) J + — In xi + — In x 2 . 

The vortex moving in the first quadrant has the intensity V and coordinates 
(x\; 7/1) and the vortex in the fourth quadrant has the intensity ~r and coordi- 
nates (x2~, 7/2). The vortices located in the left half-plane have the parameters 
-r, (-xi; 7/1) and T, (-x 2 ; y%) accordingly. 

From the first integral H = E = const one can get a relationship between 
the asymptotic values y^, xioo, X2oo- the half of the distance (t/oo) between 
the elements of the vortex pairs coming from infinity along the axis Ox, and 
the same for the pairs going to infinity along the axis Oy (xioo and X200) 

£100^200 = vlo- (10) 

The dipole power occurs implicitly through the difference of asymptotics x loo 
and X2oo- Really, as it is easily seen from the fig.3, the interaction between 
the vortices gives the same distortion in the x direction for both the 1-st and 
2-nd vortices. The difference only arises due to the influence of the flow. 

In the case when there is also a source (sink) with intensity Co in the 
origin, the relationship (10) becomes 

£1002:200 = y^exp < — \ . (11) 

Note that the relation (TIT!) remains valid if there are also the multipoles of 
the higher orders. It has been taken into account that the initial asymptotics 
corresponds to the symmetric pair of vortices. In terms of the ring vortices 
this corresponds to the dipole toroidal vortex of "infinite" radius, which is 
compressed by the interaction of the ring components. Dipole flow introduces 
an asymmetry in the motion, and different velocity correspond to vortex 
rings of different radius thrown in opposite directions. Thus, within the 
framework of the dipole-toroidal model one can naturally explain both the 
very appearance of the ejections accelerated by the accretion flow and the 
observed asymmetry of the ejections in AGNs. 
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Figure 3: The scheme of motion of four vortices in the dipole flow which 
destroys the symmetry relatively of the axis Ox. It is shown the asymmetry 
of the movement of vortices in such a stream. 

Conclusion 

The model of active galactic nuclei proposed in [10] gives the possibility 
within the framework of hydrodynamics of an ideal incompressible fluid to 
study the dynamic behavior of toroidal structures and the influence on their 
movement of accretion-wind flows. 

In this paper we considered a simplified plane model, for which the ana- 
logue of a vortex ring is a pair of point vortices, which axis coincides with the 
axis of the ring. The dipole- vortex structure of the torus in the 2D model is 
represented by the two pairs of vortices with a common axes and the angular 
momenta of the opposite signs. It takes into account that the symmetric 
pairs of vortices correspond to the initial asymptotics. 

In the terms of the ring vortices that corresponds to initial dipole toroidal 
vortex of " infinite" radius which is compressed due to interaction of the ring 
components. In the absence of the background flow this problem resembles 
the classical problem of the Helmholtz vortex ring interaction with the wall 
parallel to the plane in which the vortex ring lies. The wall can be replaced 
by a mirror image of the vortex, and the problem can be reduced to the inter- 
action of oppositely rotating vortex rings. However, in our case the direction 
of rotation is opposite to the corresponding to approaching the vortex to the 
wall, which was considered by Helmholtz. (Our choice corresponds to the 
moving off the vortex from the wall.) In a purely radial flow it admits a 
Hamiltonian formulation and exact solution of the dynamic problem [11, 12]. 
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Figure 4: Scheme of the motion of ring vortices in a converging flow, ap- 
propriate to the discussed planar analog (symmetrical flow [12]). 
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In this study the conclusion about acceleration of the ejections by the radial 
accretion flow is expanded on the general case of two-dimensional flow. It is 
shown that the monopole component of the flow is only responsible for the 
acceleration of ejections, and a dipole flow component may be responsible for 
the asymmetry of bilateral ejections. 

Thus, the dipole-toroidal model can naturally explain both the very emer- 
gence of ejections, accelerated by the accretion flow in the case of the general 
character of the flow, and the observed asymmetry of the ejections of active 
galactic nuclei and quasars. 

The text basically corresponds to [17]. We have added some figures and 
corrected typos in formulas. 

For an extended discussion of this subject see [18]. 

The authors are sincerely grateful to N.N. Kizilova for the useful notes. 
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